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We study the influence of tlie nonlinearity in tlie Schrodinger equation on the 
motion of quantum particles in a harmonic trap. In order to obtain exact analytic 
solutions, we have chosen the logarithmic nonlinearity. The unexpected result of our 
study is the existence in the presence of nonlinearity of two or even three coexisting 
Gaussian solutions. 
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I. INTRODUCTION 



The nonlinear Schrodinger equation with the logarithmic nonlinearity (we use the units 

h = 1 and m = 1) 



i dMr,t) = (^~A + V{r,t) -b \ogmr,t)\ya')) ij{r,t) 



(1) 



was introduced 



Bialvnicki-Birula and Mvcielskil | long time ago to seek possible departures 



of quantum mechanics from the linear regime. The parameter b measures the strength of the 
nonlinear interaction (positive b means attraction) and a is needed to make the argument of 
the logarithm dimensionless — it plays no significant role since the change of a results only 
in an additive constant to the potential. In what follows, we shall absorb the parameter a 
into the wave function that amounts effectively to putting a = 1. 



It has bee n proven in beautiful experiments with neutron beams 



G abler et al 



Shimonvi . 



ShuU et al 



that the nonlinear effects in quantum mechanics, if they exist at all, are 
extremely small. The upper limit for the constant b was determined to be 3.3 10~^^ eV. 
Thus, the applicability of the logarithmic Schrodinger equation to the time evolution of 
wave functions seems to have been ruled out. Nevertheless this equation, owing to its 
unique mathematical properties, has been used in many branches of physics to model the 



2 



nonlinea r behavior of various phe nomena. It has bee n appli ed in the s tudv of dissipative 



svstems 



Rulian et a) 



Hernandez and Remaud| . in nuclear phvsics 



and even in geophysics 



De Martino et al 



Hefteil . in optics 



Krolikowski et al 



In contrast to the properties of 



other nonlinear equations, the logarithmic Schro dinger equation in anv number of dimensions 



Bialvnicki-Birula and Mvcielskil ] . Gaussons 



possesses analytic solutions, called Gaussons in 
represent localized nonspreading solutions of the Gaussian shape. The internal structure of 
the Gaussons may also change in time. The existence of these analytic solutions enables one 
to study in detail the influence of nonlinearities. In this paper we focus our attention on 
the behavior of the solutions of the logarithmic Schrodinger equation in a rotating harmonic 
trap. The aim of our study was to see to what extent the nonlinear interaction may change 
the dynamics and affect the stability of solutions. Perhaps, our results will help to better 
understand the behavior of the B ose- Einstein c onde nsate in a rotat ing trap. Previous studies 



of these problems (for example. 



Recati et al, 



and 



Cozzini et al.l |) were often based on the 



hydrodynamic equations and we plan in the future to express our results in terms of the 
hydrodynamic variables. 

II. FORMULATION OF THE PROBLEM 



The logarithmic Schrodinger equation in a rotating trap has the form 
tdMr,t) = ('-iA + ir-y(t)-r-Mog(|^(r,t)|2)')^(r,t), 



(2) 



where the symmetric 3x3 matrix V{t) depends on time due to rotation. In order to simplify 
the analysis of stability, we assume that the trap is subjected to a uniform rotation and we 
shall use the coordinate system co-rotating with the trap. In this manner the potential 
becomes time-independent but due to rotation there appears an additional term in the 
equation. 



idtip{r,t) 



-Ia + ir-r-r - b\ogMT,t)\') - n-M ] z^(r,t). 



(3) 



where f2 is the vector of angular velocity and M = r x p is the operator of angular momen- 
tum. We shall seek the solutions of Eq. (jH)) in the Gaussian form 



V^(r, t) = A^(t)e^-^W exp ( -If . + iB{t))-i{t) + iTz{t)-r 



1 



(4) 
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where f = r — ^(t). The time-dependent vectors ^(t) and 7r(t) specify the position and 
momentum of the center of mass of the Gaussian wave packet and the time-dependent real 
symmetric matrices A{t) and B(t) specify the shape and the internal motion of the wave 
packet, respectively. The two real functions N{t) and f{t) define the normalization and 
the overall phase of the Gausson. Substituting this Ansatz into Eq. (jHl), we arrive at the 
following set of ordinary differential equations for all the functions entering our formula (jU 

dA{t) 



dt 
dI3{t) 
dt 

dm 

dt 
dn{t) 

dt 

dN{t) 

dt 

dm 

dt 



n,A{t) 



B{t)A{t) + A{t)B{t) 
B{t)^ - A{tf + V + 2bA{t) - tl,B{t) 

7rit)-nxm), 



-TT{B{t)}N{t), 



--[Tr{A{t)} + 7r{t)Mt)~m-V-m 



(5) 
(6) 
(7) 
(8) 

(9) 
(10) 



where the antisymmetric matrix Cl and the components of the angular velocity vector ft are 
related through the formula Qij = eijkft^ . Note, that the internal motion (described by A{t) 
and B{t)) completely decouples from the motion of the center of mass (described by ^(t) and 
7r(t)). In turn, the equations for the normalization factor and the phase can be integrated af- 
ter the internal and the c enter of mass motion has been determined. This decoupling follows 
from the general theorem 



Garcia-RipoU et al. 



and 



Bialvnicki-Birula and Bialvnicka-Birula 



stating that from every solution of a nonlinear Schrodinger equation in a harmonic poten- 
tial (including time-dependent potential) one may obtain a solution displaced by a classical 
trajectory fully preserving the shape of the wave function. 



III. SOLUTIONS AND THEIR STABILITY 



In what follows, for simplicity, we shall assume that the trap rotates along one of its 
principal axis. In this case the motion in the direction perpendicular to the rotation plane 
decouples and we are left with a two-dimensional problem. In the stationary state of our 
system the center of mass motion must be absent (^(t) = 0, 7r(t) = 0) The stationary state 
of the system is described by the wave function characterized by the solution of the following 
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two time-independent equations for two 2x2 matrices A and B 



= BA + AB 



Q,A 



= B^-A^ + V + 2bA 



Q,B 



(11) 
(12) 



We shall seek the solutions of these equations in the coordinate frame in which the matrix 
V is diagonal, V = Diagjo;^, cul}. We assume, for definitness, that ui < uji- It follows 
from Eqs. (fTTHT^ that in this frame the matrix A is also diagonal and the matrix B is 
off-diagonal. Finally, we are left with three equations for two matrix elements ai, of A 
and one matrix element (3 of B 

(«! + a2)/3 - («! - ^2)^ = 0, (13) 
0^ -a\^u)\^ 2bai + 2(3VL = 0, (14) 
p'^-al + Lul + 2ba2 - 2(3VL = 0. (15) 

It follows from Eq. ()13p that in the absence of rotation (3 must vanish and we obtain im- 
mediately two physically acceptable solutions of the decoupled quadratic equations for the 
parameters a 



a2 



+ by ujj + b), (16) 

iu2^1+ by Lui + b), (17) 

f3 = 0. (18) 

The two remaining solutions yield negative values of the a's and must be rejected. Thus, 
in the absence of rotation the nonlinearity modifies only the size of the Gaussian wave 
function without introducing any significant changes. Even for negative values of b (nonlinear 
repulsion), stable solutions described by (fT^ and (fTTj) always exist, no matter how strong 
is the repulsion. 

Simple analytic formulas can also be obtained in the presence of rotation but without 
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nonlinearity. The formulas for the Gausson parameters read in this case 



V 1 ' z ' V \ i I \— Z I , . 

"1 = . . n r... u r... ' (19) 









1 + ^{ujI - n^)/{uji 












1 + - 







«2 = . . , . , (20) 



^ . ^^-^(^iZfglMzS. (21) 
1 + v/(^i-^')/(^?-^') 

The values of a are real in the two regions of stability when Vt < uoi (region 1) and 
> UJ2 (region 2). The same regions of stability were obtained in the analysis of 
he characteristic frequencies in classic al or quantum-mechanical center-of-mass motion 



Bialvnicki-Birula and Bialvnicka-Birulal | . In the formulas (fT^ and (f^UI) the + and - sign is 



to be chosen for the region 1 and the region 2, respectively. 

In the presence of both rotation and nonlinearity the properties of solutions change sig- 
nificantly. The most striking difference is the appearance of additional stationary Gaussian 
solutions. This is an unexpected result because in the linear theory a purely Gaussian shape 
always is found for only one fundamental state of the system — all other states have poly- 
nomial prefactors. We have not been able to find closed expressions for the parameters 
a and /5, so we had to resort to numerical analysis of the solutions of Eqs. p3flT3j) . We 
present our results in three plots showing the calculated values of the parameters ai and 
a2 that determine the shape of the Gaussian wave function. These values are plotted as 
functions of the angular velocity Vt. In all plots we have fixed the trap parameters to be 
LOi = a/2/3, LO2 = a/4/3. We have chosen three values of h to describe the following charac- 
teristic cases. In Fig. ^ we plot the values of a's without the nonlinear interaction (6 = 0). 
In Fig. 121 we added the attractive nonlinear interaction (6 = 1) and in Fig. El the repulsive 
nonlinear interaction (6 = — 1). 



IV. CONCLUSIONS 



Knowing the exact analytic form of the solutions of our nonlinear Schrodinger equation we 
were able determine the influence of rotation and nonlinearity on the stability of solutions. 
The unexpected result of our analysis is that the repulsive interaction expands the region of 
stability. We have to admit, however, that this may be true only for the special form of the 
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nonlinearity: the logarithmic nonhnearity. 
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FIG. 1: This plot shows the values of ai and 02 in the absence of the nonlinearity. For each value 
of 0, in the stability regions there is just one Gaussian wave function whose shape is described by 
the values of a's. 




0,0 0,4 0,8 1,2 1,6 2,0 2,4 2,8 3,2 3,6 



FIG. 2: This plot shows the values of ai and 02 in the case of attractive nonlinear interaction. 
For sufficiently large values of b, as in this case, there are no regions of instability. For small and 
for large values of there is just one Gaussian wave function but for intermediate values there are 
two or even three solutions. Matching pairs of a's are distinguished by the lines of the same style: 
solid, dashed, and dotted. 
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FIG. 3: This plot shows the values of ai and 02 in the case of repulsive nonlinear interaction. 
The upper region of stability is extended now downwards as compared to the case without the 
nonlinear term. Moreover, there are two solutions (solid and dashed lines) that coexist in the 
newly established region of stability. 
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